Double charm production e~^e' 



> J /ip -\- cc at B factories with next-to- leading order 
QCD corrections 
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The inclusive J/ip production in e'^e~ — > J/ipcc at B factories is one of the most challenging 
open problems in heavy quarkonium physics. The observed cross section of this double charm 
production process is larger than existing leading order (LO) QCD predictions by a factor of 5. 
In the nonrelativistic QCD (NRQCD) factorization formalism, we calculate the next-to-leading 
order (NLO) QCD virtual and real corrections to this process, and find that these corrections can 
substantially enhance the cross section with a K factor of about 1.8. We further take into account 
the feeddown contributions from higher charmonium states (mainly the i^{2S) as well as Xc./) and 
the two-photon contributions, and find that the discrepancy between theory and experiment can be 
largely removed. 
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The double charm production in e^e~ annihilation at 
B factories is one of the most challenging open problems 
in heavy quarkonium physics and nonrelativistic QCD 
(NRQCD) (for a review of related problems, see[l|). The 
exclusive production cross section of double charmonium 
in e+e^ J/i^Vc at = 10.6 GeV measured by Belle 
[3,[| and BaBar[| are larger than the leading order (LO) 
calculations 3.8 ~ 5.5fb [1, i] in NRQCD by possibly 
almost an order of magnitude (see also [7[). 

Moreover, the inclusive J/i(j production cross section 



via double cc in e"' 
sured by Belle 



J/^jjcc at = 10.6 GeV mea- 



a[e+e- J/iP + cc + X] = (CSTl^'ig ± 0.17) pb , (1) 

is about a factor of 5 higher thanjjheoretical predictions 
including both the color-singlet ^ |9| and color-octet 9] cc 
contributions at leading order(LO) of as in the NRQCD 
factorization formalism [IO5. This is another intriguing 
challenge in the double charm production problem, aside 
from the exclusive J/tpVc production in e'^e^ annihila- 
tion. 

Some theoretical studies were attempted in order to re- 
solve the large discrepancy in e"*" -I- e~ — *■ J/ip + cc. Liu, 
He, and Chao calculated the color-octet contribution Q 
and the two-photon contribution to J/ip + cc produc- 
tion 111 in NRQCD. But those contributions are small 



and can not make up such a large discrepancy. Hagiwara 
et al. assumed a large renormalization K factor {K ~ 4) 
for the J/ipcc cross sectionfl2i]. Kaidalov introduced a 
nonperturbative quark-gluon-string model 13]. Kang et 
al. got cr(e+e" J /ip + cc + X) / a{e+ J/ip + X) = 
0.049 in the color-evaporation-model Other sugges- 
tions to resolve this problem may be found in Ref. [l|. 

In order to further clarify this problem, in this paper 
we present a result for the next-to-leading order (NLO) 
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FIG. 1: Two of the four Born diagrams for e~ 
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QCD correction to the inclusive J/-0 production pro- 
cess of 6+ -I- e~ J/i^ + cc. And we have already 
found that the NLO QCD correction to the exclusive 
production process of e+ -I- ~> J/i^ + Vc is crucial, 
for which the K factor (the ratio of NLO to LO ) may 
reach to a value of about 2, and hence essentially reduces 
the large discrepanc y b etween theory and experiment of 
e+ -I- e~ J/^ + ricuM (with significant relativistic cor- 
rections further considered this discrepancy is probably 
resolved [lii [l7|: enhancement effects due to use of light- 
cone formalism and relativistic corrections are also pro- 
posed in [lit). 

At LO in as, J/ip + cc can be produced at or- 
der a'^a'j. (see, e.g., Refs. [1, [11]). There are four 
Feynman diagrams, two of which are shown in Fig. [l] 
and the other two diagrams can be found through re- 
versing the arrows on the quark lines. Momenta for 
the involved particles are assigned as e^(fci)e+(fc2) 
1*{Q) ~^ J/i'{pj) + c(pc) + c{pc). In the calculation, 
we use FeynArts [l9| to generate Feynman diagrams 
and amplitudes; FeynCalc (20| for the tensor reduc- 
tion; and LoopTools [2H for the numerical evaluation 
of the infrared (IR)-safe one-loop integrals. Finally we 
use Mathematica to integrate over phase space and get 
the numerical results. 

To NLO calculation, the cross section is 

= CTBorn + Cr virtual + (y real + 0{a^a1.), (2) 
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where 



da 



virtual 



i ^ E 2 Re{M*B,„,MNLo)dPS3, 

dareal = j ^ l"^-"' I'^P^^- (3) 

Here the factor l/2s is the flux factor. ^ means sum 
over the polarizations of initial and final state particles. 
dPSa and dPS4 are the three and four-body phase spaces 
respectively. 

There are ultraviolet (UV), infrared(IR), and Coulomb 
singularities, and we treat them in the same way as 
in Ref. 15|. For the box diagrams shown in Fig. [21 



BoxN8 and BoxNlO have IR and Coulomb singularities, 
BoxN3 does not have IR singularity, while the other 
nine diagrams have IR singularity. BoxNl + BoxN4 
, BoxN6 + BoxN7 + PentagonN12 are IR finite respec- 
tively. IR terms of BoxN9 -I- BoxN2 -t- PentagonNll are 
canceled by Vertex diagrams. The IR term in counter 
terms and BoxN5 -I- BoxN8 -I- PentagonNlO should be 
canceled by the real corrections. And the Coulomb singu- 
larity terms in BoxN8 + PentagonNlO should be map ped 
into the wave functions of J/V'- Similar to Ref. [15[, we 
use Z? = 4 — 2e space-time dimension and the relative ve- 
locity V to regularize the IR and Coulomb singularities. 
The IR and Coulomb singularity terms in the virtual cor- 
rections are 



da. 



I R. Coulomb 
virtual 

Btt I 



1 _ Pc ■Peln(-Seg) 1 



-(4) 



•y/l-4m2/(pc-|-pc)^ 



where Sr, , 

^l-4m2/(p,-HpE)2 + l 

There are thirty diagrams for real corrections, and half 
of them are shown in FiglS] The other fifteen diagrams 
can be obtained through reversing the arrows on the 
quark lines that are connected with J 1^. The calculation 
of the real corrections is similar to the leading order cal- 
culation, but there should appear the IR singularity [i^] • 
We find that Real Nl, Real N7, Real N12, and Real N14 
are associated with IR singularity. And the eikonal fac- 
tors of Real N4, Real N6, Real N8, Real Nil, Real N13, 
and Real N15, in which the gluon is connected with the 
external charm quark and anti-charm quark in the J/V', 
are canceled by themselves. The other five diagrams 
Real N2, Real N3, Real N5, Real N9, and Real NIO are 
independent of IR singularity. Using the eikonal approxi- 
mation, we get the IR singularity terms in real corrections 



IR 

real 



\n{-Scc)Pc ■ Pc 



(5) 



The Coulomb singular terms in c'o'^^j^""'"™'' can be 
mapped into the J/ip wave function. 

We now turn into numerical calculations for the cross 
section of e+ -I- e~ J/tp + cc + X. To be consistent 
with the above NLO calculation, the value of the J/ip 
wave function squared at the origin should be extracted 
from the leptonic width at NLO of as (see e.g. (lo| ) 



16^2 (1 - 4:Cfois/-k) 



T{J/xP e+e-). (6) 



Using the experimental value 5.55±0.14±0.02KeV [23|, 
we obtain |i?s(0)p = l.OlGeV^, which is a factor of 1.25 
larger than 0.810GeV^ that was used in Ref. 9] from po- 
tential model calculations. Taking mj/^ — 2m (in the 

nonrelativistic limit), m = 1.5 GeV, A^-l^ = 338MeV, 

^' 'MS ' 

cks(/^) = 0.259 for fj, — 2m (these are the same as in 
Ref. 9] except here a larger |i?5(0)p is used), the cross 
section for e+ + J/il) + cc + X at NLO of as is 



(7(e^ 



J/ip + cc + X) = 0.33 pb. 



(7) 



It is a factor of 1.8 larger than the LO result 0.18 pb 
obtained with the same parameters. For fi — m and 
\/s/2, we have as = 0.369 and 0.211, and get the cross 
section 0.53 pb and 0.24 pb respectively. If we set m = 
1.4GeV and /i = 2m, the cross section at NLO of as is 



a(e^ 



J/ip + cc + X) = 0.47 pb. 



(8) 



It is about a factor of 1.7 larger than leading order cross 
section 0.27 pb. The dependence of the cross section on 
the renormalization scale is shown in Fig. U) When n 
changes from mc = 1.5 GeV to \/s/2 — 5.3 GeV, the 
ratio a^jj) / a{y/s /2) is found to vary from 3.05 to 1 in 
LO result. For NLO, with TOc = 1.4(1.6) GeV, the ratio 
a{n)/a{y/s/2) varies from 2.22(2.26) to 1. We see that, 
as expected, the scale-dependence in NLO is consider- 
ably reduced compared with that in the LO. A detailed 
discussion will be given elsewhere. 

We should also include the QED contribution as well 
as the two-photon contribution of e+e~ 27* >^/V' + 
cc. Furthermore, since the experimental data are for the 
prompt J -\- CC + X production, we should consider the 
feeddown contributions from higher charmonium states 
such as e+e^ — > ip{2S) + cc + X J/4> + cc + X, and 
e+e" XcJ + cc ^ J/ip + CC + X. 

Two of the six QED diagrams of e+e^ ^ j* ^ J/'4' + 
cc are shown in Fig. [5l The other four diagrams can be 
obtained by replacing the gluon with the photon shown in 
Fig. [TJ Contributions from QED diagrams can interfere 
with that from QCD Born diagrams, and resulted in a 
cross section of 8 fb at order 0{asa^). 

e+e~ 27* J/ipcc has been calculated by Liu 
et al. 11|. Using their result, this cross section is 23 x 



They just cancel the IR singular terms of da^^'^"]^''""^^ . 



fb = 29 fb, where the factor ^j^o is due to using 
the new value of the J/iJj wave function at the origin. 
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FIG. 2: Twelve of the twenty-four box and pentagon diagrams for e~(fci)e^(fc2) J /'4'{pj) + c(pc)c(ps)- 




Real Nil Real N12 Real N13 Real N14 

FIG. 3: Fifteen of the thirty real correction diagrams for e~ J/i/' + cc+ ( 
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At leading order in v (the relative velocity of quark 
and antiquark in the charmonium rest frame), the differ- 



ence between e+e V'(25') -|- cc and e+e 



cc 

is in the wave functions at the origin. Using Eq.®, the 
contribution from the transition ip{2S) J/tp is to en- 
large the direct production of J/tp -I- cc + X by a fac- 
tor of ^2|g}pB(^(25) ^ J/^X). By using Ti^{2S) ^ 
e+e^) = 2.48 ± 0.06KeV and the branching ratio for the 
ip{2S) J/ipX transition fraction B = 56.1 ±0.9% [H], 
we find the enlarging factor to be 0.355. 

The cross sections of e+e^ — * XcJ + cc + X with 
both color singlet and octet contributions were calcu- 
lated in Ref.[9|]. Moreover, the color octet contribu- 
tion to J/ipcc production was also estimated to be 



about 11 fb in Rcf. 9]. Using their results and the 
observed branching ratios for Xcj ^ J li^l transitions 
B = 1.31%, 35.6%, 20.2% for J = 0, 1, 2 respectively H, 
we find the sum of feeddown from Xcj and color-octet for 
J /if: contributions to be 21 fb. 

Combining all these contributions, the prompt cross 
section of e+e" — > J/tj} + cc + X at NLO of as is 



(7 



prompt 



X) = 0.51pb. (9) 



It is 59% of the data value in Eq. ([T]). If we set m — 
1.4GeV and /i — 2m, and ignore the differences of other 
contributions due to the change of mass, then the prompt 
cross section of e+e^ J/tp + cc + X at NLO of as is 

<Jpromptie+ + e- ^ J/i/j + cc + X)^ 0.71 pb. (10) 
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FIG. 4: Direct production cross sections of e"'" + e~ 
J/ijj + cc+X as functions of the renormalization scale /i. Here 
|7?s(0)p = 1.01GeV^ A = 0.338GeV, ^ = 10.6GeV; NLO 
results are represented by solid lines and LO one by dashed 
lines; the upper line is for m = 1.4GeV and the corresponding 
lower line is for m = l.SGeV. 




FIG. 5: Two of the six QED diagrams for e'^e — > J/iIjcc. 



It is 82% of the experimental value in Eq. ([T]). 

In order to see the uncertainties due to parameters 
TOc, fi, as{n), we set fi — 2.8(5.3) GeV, m = 1.4 =f 
0.2 GeV, and then get the cross section in the range 
0.71 (0.53;^:^^) pb. The NLO QCD correction to 
e+ + e~ J/ip + cc is large, despite of existing paramet- 
ric uncertainties. 

In conclusion, we find that by taking into consideration 
all NLO virtual corrections with self energy, triangle, box, 
and pentagon diagrams, and the real corrections, and fac- 
toring the Coulomb singular term into the cc bound state 
wave function, we get an ultraviolet (UV) and infrared 
(IR) finite correction to the direct production cross sec- 
tion of e+e" — > J/ip + cc at -^/i = 10.6 GeV, and find 
that the NLO QCD correction can substantially enhance 
the cross section with a K factor of about 1.8. With 
TO = 1.4GeV and fi — 2m, the cross section of direct 
J/ipcc production through one-photon is estimated to be 
0.47 pb. Adding the feeddown contributions from higher 
charmonium states (mainly the V'(2S') as well as Xcj) 
and contributions from two-photon process and color- 
octet channels, the prompt production cross section of 
e+e~ J/ tp + cc at NLO in is found to be 0.71 pb, 
which is 82% of the experimental value 0.87 pb. Hence 
the discrepancy between theory and experiment is largely 
removed, despite of certain theoretical uncertainties. 
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